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Abstract. We study global controllability of 'rotating' multidimen- 
sional rigid body (MRB) controlled by application of few torques. Study 
by methods of geometric control requires analysis of algebraic structure 
introduced by the quadratic term of Euler-Frahm equation. We discuss 
problems, which arise in the course of this analysis, and establish several 
global controllability criteria for damped and non damped cases. 



1. Introduction 

In recent work [H [U [3] one studied controllability of Navier-Stokes (NS) 
equation, controlled by forcing applied to few modes on a 2D domain. Geo- 
metric control approach has been employed for establishing approximate 
controllability criteria for NS/Euler equation on 2D torus, sphere, hemi- 
sphere, rectangle and generic Riemannian surface with boundary. 

In the present contribution we address controllability issues for a finite- 
dimensional "kin" of NS equation - Euler-Frahm equation for rotation of 
multidimensional rigid body (MRB) subject to few controlling torques 
and to possible damping. The equation evolves on so(n). We formulate 
global controllability criteria which are structurally stable with respect to 
the choice of 'directions' of controlled torques. 

According to geometric approach to studying controllability, one starts 
with a system controlled by low-dimensional input and proceeds with a 
sequence of Lie extensions ([SI [6]) which add to the system new controlled 
vector fields. The latter are calculated via iterated Lie-Poisson brackets of 
the controlled vector fields and the drift (zero control vector field). The core 
of the method and the main difficulty is in finding proper Lie extensions and 
in tracing results of their implementation. 

The Lie extension employed in [U El [3] for studying controllability of NS 
equation, and similar one used equation below (see Subsection 14. ip . involves 
double Lie bracket of drift vector field with a couple of constant controlled 
vector fields (they are identified with their values, or directions belonging to 
so(n)). At least one of the directions must be a steady state of MRB, i.e. 
an 'equilibrium points' of Euler-Frahm equation. 
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The double Lie bracket results in constant controlled vector field (ex- 
tending direction); the correspondence between couple of original controlled 
directions and the extending one defines bilinear operator (3 on so{n). More 
extending controlled directions are obtained by iterated application of j3. For 
proving global controllability of MRB we must verify saturating property - 
coincidence of the set of extending directions with so(n) after a number of 
iterations. 

Tracing the iterations is by no means easy. For NS equation all cases, 
successfully analyzed in [IJ El [3], are related to an explicit description of 
the basis of steady states and to specific representation of the operator (3 
with respect to this basis. The results, so obtained, are heavily dependent 
on choice of original controlled directions and on geometry of the domain 
where the NS equation evolves. 

Below we manage to establish several controllability criteria for damped 
and non damped MRB controlled by one, two or three torques. We pay spe- 
cial attention to deriving criteria which are structurally stable with respect 
to perturbation of (some of) the controlled directions. 

2. Euler Equation for Generalized Rigid Body and 
Euler-Frahm equation for MRB 

We follow [1] for definition of 'generalized rigid body'. Let Q be a Lie 
group, g its Lie algebra and let left-invariant Riemannian metric on Q be 
defined by scalar product (•, •) on g. 

Introduce X : g i— > g* - a symmetric operator, which corresponds to the 
Riemannian metrics by formula: (£,T)) = Z£\t), where - 1 - is the natural pair- 
ing between g and g*. The operator Z is called inertia operator of generalized 
rigid body. 

The trajectory of the motion of generalized rigid body is a curve g{t) 6 Q. 
Angular velocity, corresponding to this motion is: £1 = L g -i*g 6 g, where L g 
is left translation by g. The image of angular velocity £1 under I is angular 
momentum M G g*. Energy of the body equals (O, O) = M|0. 

Euler equation for the motion of generalized rigid body is 0, = S(fi,Q), 
where bilinear operator B : g x g i— > g is defined by formula: 

(1) ([a,b],c} = (B(c,a),b), 

[•, •] staying for Lie-Poisson bracket in g. 

MRB is particular case of generalized rigid body, where the Lie group 
Q = SO(n), and angular velocities G g = so(n) are skew-symmetric 
matrices. 

Identifying so(n) with so*(n) by means of Killing form, we may think of 

momentum M as of skew-symmetric matrix. Then the inertia operator is a 

K 

map Tc ■ !— > {QC + CO) = M € so(n) = so*(n), where C is some positive 
semidefinite matrix. 
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Operator Xc is symmetric with respect to Killing form and is invertible 
(Sylvester theorem), whenever C is positive definite. 

We compute B according to ([T]) ([•,•] being matrix commutator): 

B{Q},£i 2 ) =z c 1 [icti 1 ^ 2 ] , g(n,n) =i c 1 [i c n,n] = i^[c,n 2 ]. 

Euler-Frahm equation for the motion of MRB is: 

(2) n = i c 1 [i c n,n}=i c 1 [c,n 2 ], 

The motion, subject to damping, is described by the equation 

n = i c 1 [c,n 2 } -vn, v> o. 

3. Controllability of rotating MRB: problem setting and main 

RESULTS 

Controlled rotation of MRB is described by equation 

r 

(3) n=l c 1 [C,n 2 ]-un + ^2G i u i (t), u>0, C?€ao(n). 

1=1 

We are interested in global controllability of ([3]), meaning that for any 
f2, £l £ so(n) system ([3]) can be steered from £1 to Cl in some time T > 0. We 
are interested in achieving global controllability by small number of controls; 
we prove that r can be taken < 3 for all n > 3. 

Equation ([3]) is particular case of control-affine system with quadratic (+linear) 
drift vector field and constant controlled vector fields. 

The following genericity condition is assumed to hold furtheron: symmet- 
ric matrix C is positive definite and has distinct eigenvalues. 

Our first result claims global controllability of MRB by means of two 
controlled torques. 

Theorem 3.1. There exists a pair of directions G^G 2 6 so(n) (depending 
on C ), such that the system (G|) with r = 2 is globally controllable. □ 

The proof of this Theorem, sketched below, is based on direct computa- 
tion of Lie extensions in specially selected basis, related to C . More difficult 
is formulating criteria, which are structurally stable with respect to pertur- 
bation of controlled directions. 

We start with non damped MRB, controlled by one torque. In this case - 
given recurrence of Euler-Frahm dynamics (|2|) - bracket generating property 
suffices for guaranteeing global controllability. This property means that 
evaluations (at each point) of iterated Lie brackets of drift and controlled 
vector fields span so(n). Given high dimension of so(n), verification of 
the bracket generating property for generic controlled direction is nontrivial 
task. We do this analyzing linearization of quadratic Euler operator. The 
result is 

Theorem 3.2. For generic G € so(n) the system Cl = 1^}[C, Q 2 ] + Gu(t), 
is globally controllable, also if control is bounded: \u\ < b, b > 0. □ 
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We now pass to the damped case. Our method requires one of the con- 
trolled directions to be steady state for MRB. Recall that steady state 
or steady direction of MRB is equilibrium point of ([2j) - a matrix G for 
which [2 C G, G] = [C, G 2 ] = 0. Matrix G is principal axis of MRB, if 
XqG = fiG, fj, £ R. These two sets coincide for n = 3, while for n > 4 the 
set of steady directions is much richer. 

The results obtained for the damped case differ for odd and even n. 

Theorem 3.3. Let r = 2, n be odd in For generic stationary direction 
G l and generic G 2 € so(n) the system {3|) is globally controllable. □ 

An additional symmetry in the case of even n, obliges one to involve 
additional controlled direction for achieving global controllability. 

Theorem 3.4. Let r = 3, n be even in (E|j. For generic stationary direction 
G 1 £ so(n) and generic pair (G 2 , G 3 ) of directions the system is globally 
controllable. □ 

Generic element of a subset W C so{n) means an element of open dense 
subset of W in induced topology. 

4. Sketch of the proof of Theorem 13.11 

4.1. Key Lie extension. Lie extensions mean finding vector fields X, 
which are compatible with control system, in the sense that closures of at- 
tainable sets of the control system are invariant for X. If one is able to 
prove global controllability of the system extended by some compatible vec- 
tor fields, then controllability of the original system can be concluded by 
standard argument. 

Key Lie extension, we employ, is described by the following 

Proposition 4.1. Let for control system 

(4) x = f(x) +g(x)u + g(x)v, 
evolving on a manifold Q, hold the relations 

(5) {~g,g} = 0, {g,{9,f}} = 0, 

({■,■} stays for Lie brackets of vector fields on Q). Then the system x = 
f(x) + g(x)u + g(x)v + {g, {g, f}}(x)w is Lie extension of □ 

Remark 4.1. Vector fields ±{g,{g,f}} are extending controlled vector 
fields; they are also compatible with Q). □ 

We will repeatedly employ Proposition 14.11 for extending control system 
(|3|). At each step the first of the relations ([5]) will be trivially satisfied 
since all original and extending controlled vector fields will be constant. For 
drift vector field f(£l) = I^, l [C, Q 2 ] in ([3]), and constant controlled vector 
field g = G G so(n) , the Lie bracket {g, {g,f}} = T~^}\C, G 2 ] is constant 
vector field. The second relation ([5]) would hold if and only if G is steady 
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state. When repeating the extension it is crucial to guarantee at each step 
disponibility of steady state controlled direction. 

For two constant controlled vector fields g = G,g = G, G, G € so(n) the 
value of constant extending controlled vector field {g, {g,f}} is 

(6) (3(G,G) =1~ 1 [C,GG + GG}; 

([6|) defines symmetric bilinear operator (3 on so(n). 

4.2. Algebra of principal axes and controllability proof. Diagonalize 
matrix C presenting it as C = Ad SD = SDS' 1 with S orthogonal and 
D = diag{7i, ...,/„}, h < I 2 < ■ ■ ■ < I n - 

Introduce matrices @ rs = l rs — l sr € so(n), (1 < r < s < n), with l rs 
being matrix with (the only nonvanishing) unit element at row r and column 
s. Matrices £l rs = Ad SQ rs turn out to be 'eigenvectors' of the operators 
(ad C) and Ic- They form set of principal axes of the MRB. 

'Multiplication table' for (3 with respect to the basis VL rs is 

P(n rs , n rs ) = o, 0(n rs , n re ) = (i £ - i s )(i s + h)- 1 ^, 

(3(Q rs , Q ke ) = 0, whenever r, s, k,£ are distinct. 
Take 

G 1 = Q 12 - principal axis, G 2 = Q 23 + O 34 H h O"- 1 '". 

It suffices to prove that iterated applications of f3 to G 1 , G 2 result in a basis 
of so(n), because then the extended system would possess full-dimensional 
input and therefore would be globally controllable. The original system ([3]) 
would be globally controllable as well. 

According to the multiplication table G 3 = j3(G 2 ,G l ) = /3(0 12 ,0 23 ) co- 
incides up to a multiplier with principal axis Q 13 . Calculating subsequently 
extending controlled 'directions' G % = [5(G l ~ l , G 2 ), i > 2, we see that all G % 
coincide up to a nonzero multiplier with O 1 '*, i.e. are principal axes. Also 
/3(^ lj , Q lk ) coincides up to a multiplier with this means that iterating 
applications of /3 to G , G 2 generate basis of so(n). 

5. Concluding remarks 

1. As one can see proof of Theorem 13.11 is "rigid construction", based 
on specific choice of controlled directions and on computation of iterated 
Lie extensions with respect to specific basis of principal axes of MRB. If 
one perturbs one of the original controlled directions the constructions fails, 
as far as first Lie extension does not result in new stationary direction of 
MRB, and the Proposition 14 . 1 1 can not be iterated. 

This rigidity of controllability criteria with respect to the choice of con- 
trolled directions, manifested itself also in previous study of approximate 
controllability of NS system on particular 2D domains ([HE]). It does not 
seem natural, and is rather related to the proposed method. 
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Indeed one would expect structural stability of controllability criteria and 
this is achieved in the formulations of Theorems I3.2|3.3|3"74| which are struc- 
turally stable with respect to the choice of (some of) the controlled direc- 
tions. The method for establishing these criteria differs from the previous 
one. It is based on study of linearization at a steady state of Euler operator 
for MRB. The proofs will appear elsewhere. 

Besides its interest for studying controllability of MRB, the method can 
be extended onto infinite-dimensional case, and be applied to controlled 
NS/Euler equation for fluid dynamics on general 2D and 3D domains. The 
results will appear in further publications. 

2. Publication [8] studies controllability of non damped MRB by using 
of a pair of controlled 'flywheels' - different type of "internal-force controls", 
with dynamics described by bilinear control system on Lie group. 
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